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Abstract. Solar wind electrons are accelerated and reflected upstream 
by the terrestrial bow shock into a region known as the electron foreshock. 
Previously observed electron spectra at low energies are consistent with a 
fast Fermi mechanism, based on the adiabatic conservation of the magnetic 
moment (/i) of the accelerated electrons. At higher energies, suprathermal 
power law tails are observed beyond the level predicted by fast Fermi. 

The SWEA and STE electron detectors on STEREO enable measurements 
of foreshock electrons with good energy resolution and sensitivity over the 
entire foreshock beam. We investigate the electron acceleration mechanism 
by comparing observed STEREO electron spectra with predictions based on 
a Liouville mapping of upstream electrons through a shock encounter. The 
foreshock electron beam extends up to several tens of keV, energies for which 
the Larmor radii of electrons is tens of km or greater. These radii are com- 
parable to the scale sizes of the shock, and /i conservation no longer applies. 
We show that the observed enhancement in the foreshock beam beyond fast 
Fermi levels begins at energies where this assumption breaks down. 
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We also demonstrate, using the Liouville mapping technique, that the strahl 
plays an important role in the formation of the bump on tail instability. We 
discuss this in the context of recent observations in the foreshock and in so- 
lar wind magnetic holes. 
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1. Introduction 

The region of interplanetary space which is upstream of the terrestrial bow shock and 
magnetically connected to the shock is known as the foreshock. Observations of upstream 
electron beams [Anderson et al, 1979; Feldman et al, 1983; Gosling et al, 1989] near 
the leading edge of the foreshock demonstrated that incident solar wind electrons are 
accelerated by the quasiperpendicular bow shock back into the interplanetary medium. 
These electron beams are frequently associated with electrostatic Langmuir oscillations. 
Filbert and Kellogg [1979] pointed out that the effect of velocity dispersion of accelerated 
electrons convecting downstream with the solar wind is to create bump on tail electron 
velocity distribution functions (eVDFs). The bump on tail eVDFs are then unstable to 
the growth of Langmuir waves via Landau resonance. 

The canonical picture of the quasiperpendicular electron acceleration process is the 
Fast Fermi model developed by Leroy and Mangeney [1984] and Wu [1984], which treats 
the electron-shock encounter assuming conservation of the first adiabatic invariant (/i = 
mv\/2B) in the de Hoffmann- Teller (HT) frame. In the HT frame, the upstream bulk 
velocity u u and magnetic field B u are parallel, and the electrons do not encounter a 
motional electric field. The transformation velocity v ht to the HT frame from any shock 
stationary frame lies in the plane of the shock and is given by 



In this picture, the energization is a consequence of the boost to the HT frame, the 
mirroring, and the boost back to the solar wind frame. The accelerated electrons gain 



v H t — 



n x (u u x Bu) 
B u n 



(1) 



DRAFT 



February 17, 2012, 4:16am 



DRAFT 



PULUPA ET AL.: STEREO FORESHOCK ELECTRONS X - 5 

parallel velocity proportional to 2-uht- Detailed expressions for the eVDFs upstream of 
the shock in the foreshock region have been presented by Cairns [1987] and Fitzenreiter 
et al. [1990]. A review of the physics of electron acceleration at the bow shock is presented 
in Burgess [2007]. 

The fast Fermi theory works well to describe relatively low energy electron acceleration. 
Observed features of foreshock electrons such as prominent bump on tail eVDFs [Fitzenre- 
iter et al, 1984, 1996] and energy dependent loss cone widths [Larson et al, 1996] can be 
explained in the context of the theory and used to deduce information about the shock at 
the acceleration site. Simulations of electron interactions with planar and curved shocks 
[Krauss-Varban et al, 1989; Krauss-Varban and Burgess, 1991] also agree well with the 
fast Fermi model. 

However, for the more energetic electrons observed in the foreshock (for example, the 
> 16 keV electrons in Figure 2 of Anderson et al. [1979], and the power law tails observed 
by Gosling et al. [1989] and Oka et al. [2006]), the assumptions used in the fast Fermi 
model are inapplicable. Electrons with energies of tens of keV or higher have Larmor radii 
tlc which are comparable to the ion Larmor radius ru and/or the ion inertial length c/u p \ 
in the solar wind at 1 AU. Since the scale size of the quasiperpendicular bow shock is set 
by these ion scales [Bale et al, 2003], the conditions for ji conservation are violated, and 
the presence of the observed suprathermal electrons must be explained by physics beyond 
single encounter fast Fermi acceleration. 

This study presents STEREO observations of foreshock electron beams with good energy 
resolution over the entire energy range of the beams. Section 2 describes the STEREO 
instrumentation, and Section 3 describes the fitting procedure for the incident population 
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of solar wind electrons. Section 4 discusses the method of Liouville mapping the incident 
eVDF through the shock encounter. The Liouville mapping predictions are compared 
with observations in Section 5, and the impact of the results on the theory of electron 
shock acceleration are discussed. Section 6 summarizes the conclusions of this study. 

2. STEREO Measurements 
2.1. Overview 

During the early part of the STEREO mission, both spacecraft spent significant time in 
the terrestrial electron foreshock. Upstream Langmuir waves [Malaspina et ai, 2009] and 
electron beams were observed over distances of several hundred -Re upstream of the shock. 
This paper focuses on one foreshock event, which occurred at approximately 0730 UT on 
20 December 2006. At this time, the spacecraft GSE coordinates were approximately 
[77.9, -64.4, -10.8}R E . 

Langmuir wave observations were made with the S/ WAVES LFR instrument [Bougeret 
et ai, 2008], which measures electric field fluctuations from 2.5 to 160 kHz, a range which 
includes the local plasma frequency. Magnetic field measurements were made with the 
STEREO/IMPACT MAG instrument [Luhmann et ai, 2008; Acuna et ai, 2008]. Solar 
wind and foreshock electron observations were made with the STEREO/IMPACT electron 
detectors. The instruments used were the electrostatic SWEA detector [Sauvaud et ai, 
2008], which measures electrons from 1 eV to 2 keV with nearly 4ir angular coverage, and 
the solid state detector STE [Lin et ai, 2008], which measures electrons from 2 keV to 
100 keV. Together, these detectors span the energy range of the known populations of 
quiet time solar wind electrons [Lin, 1998]. 
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An overview plot of the foreshock event is shown in Figure 1. From top to bottom, 
Figure 1 shows the magnetic field in GSE coordinates, the LFR electric field spectrogram, 
an energy spectrogram from the STE D3 detector, STE pitch angle distributions (PADs) 
at 13.6 keV and 5.2 keV, and SWEA PADs at 1060 eV and 400 eV. For the entire event, 
the x GSE component of the magnetic field is negative, which implies that electrons with 
pitch angles close to 0° are traveling towards the bow shock (i.e., incident to the shock) 
and those with pitch angles close to 180° are traveling sunward (i.e., backstreaming from 
the shock). 

Three time intervals are indicated in Figure 1. During the first interval (indicated by 
two solid vertical lines), from 0720 to 0722 UT, the spacecraft was in the undisturbed 
solar wind, unconnected to the bow shock. The electron PADs for this time period show 
evidence of a beam-like strahl population, apparent as an enhancement in pitch angles 
close to 0° in the 400 eV PAD. 

The spacecraft then enters the electron foreshock. Starting just before 0725 UT, fore- 
shock electron beams appear as enhancements in pitch angles close to 180°. The foreshock 
electron beams are visible in the SWEA and STE PADs. Bursty beam-driven Langmuir 
wave emission at around 19 kHz is also apparent. The STE energy spectrogram shows that 
the foreshock electron beam evolves rapidly. This evolution is an expected and previously 
observed consequence of the time-of-flight velocity dispersion in the foreshock [Filbert and 
Kellogg, 1979]. 

The second and third time intervals (indicated in Figure 1 by vertical dashed lines) 
show observations at two different foreshock depths. During the first interval at 07:25, 
the spacecraft is close to the foreshock edge, and the foreshock beam energy spectrum 
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reaches well into the STE energy range. During the second interval at 07:27, the spacecraft 
is deeper in the foreshock and the foreshock beam is at lower energies. These intervals 
will be examined in detail in Section 5. 

2.2. Interpretation of STE observations 

The STE instrument has two sensor units, each consisting of a set of four thin window 
silicon semiconductor detectors. The two units are centered along the nominal Parker 
spiral in the upstream (STE-U) and downstream (STE-D) directions. [Lin et al, 2008; 
Wang et al, 2010]. 

Because the central look detection of each STE angular bin lies in the ecliptic, the 
pitch angle coverage of the STE instrument depends on the magnetic field observed at the 
spacecraft. For example, if the field is entirely in the zqse direction, each STE detector 
observes the same pitch angle (90°). If the field is in the x GS e — Ygse plane, more pitch 
angle information can be recovered. Evidence of this dependence is visible in the STE 
PAD plots in Figure 1, which demonstrate considerably wider range in pitch angle as the 
z gse component of the field decreases. 

When STEREO is in its normal orientation with the SECCHI imagers [Howard et al, 
2008] looking at the sun (which is the case for the entire mission, excepting some early 
orbital maneuvers and short, isolated magnetometer calibration intervals), multiply re- 
flected solar optical and UV light contaminates the STE-U detectors. For this reason, 
only STE-D data is presented in this paper. 

Because the detectors are open, without curved path electrostatic deflectors or foil 
covers, STE measures ions, neutrals [Wang et al, 2010], and X-rays [Hsieh et al, 2009] in 
addition to electrons. For this reason, we must be careful in interpreting STE data in the 
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foreshock region. In order to be assured that an observed STE signal is indeed a foreshock 
electron beam, we require both (a) simultaneous electron measurements from SWEA and 
(b) associated Langmuir wave activity. The STE observations themselves must also be 
consistent with electron beams. An example of STE data which is inconsistent with 
foreshock electron beams is visible in Figure 1 at 07:34 UT. Were the measured spectra 
at these times due solely to electrons, the features observable in the STE spectrogram at 
several tens of keV would be highly unstable, so the electron measurement at these times 
is likely contaminated by foreshock ions. 

The two foreshock intervals highlighted in Figure 1 have STE spectra consistent with 
electrons, simultaneous SWEA observations, and associated S/WAVES Langmuir waves. 
For the remainder of this paper the observations at these intervals are taken to be foreshock 
electron beams. 

3. Incident electrons 

The solar wind electrons comprise several distinct populations; the populations con- 
sidered in this paper are the core, the halo, the strahl, and the superhalo. The thermal 
core population makes up most of the number density of the solar wind and can be mod- 
eled well using a bi-Maxwellian distribution. The suprathermal halo population is better 
described by a bi-K distribution, which resembles a bi-Maxwellian at low energies and 
behaves like a power law distribution at high energies [Maksimovic et al, 1997]. The 
strahl component of the solar wind is a beam-like feature propagating along the magnetic 
field in the antisunward direction. A detailed description of the core, halo, and strahl 
components can be found in Stverdk et al. [2009]. 
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In addition to the above components, the quiet time solar wind electrons include a 
population known as the superhalo, which can be approximated using a power law from 
energies of several keV up to 100 keV [Lin et ai, 1996; Wang et ai, 2012]. The superhalo 
electrons have been shown to be nearly isotropic [Lin, 1998]. 

The energy spectrum and the PAD of the incident solar wind electrons are shown in 
the top row of Figure 2. Electrons from 57.5 eV to 1.1 keV are measured by the SWEA 
instrument and are marked with asterisks. SWEA exhibits a degradation in low energy 
electron measurements due to instrumental charging effects [Fedorov et ai, 2011], so the 
energies below 50 eV are not plotted in Figure 2. SWEA measurements above 50 eV are 
corrected by the transmission function described in Fedorov et al. [2011]. The incident 
PAD shown in Figure 2 shows evidence of a strong strahl beam (the peak near zero pitch 
angle) and a relatively isotropic halo population (the relatively flat lines at pitch angles 
greater than 90°). 

Electrons from 2 keV to 100 keV are measured by STE and are marked with diamonds. 
During the time interval when the incident electrons were measured, the magnetic field 
contained a strong component in the zgse direction, and the angular coverage of the STE 
instrument was limited to a narrow range of pitch angles. For clarity, not all of the STE 
electron energy bins are shown. 

Incident energy spectra at near-parallel (green) and near-perpendicular (blue) pitch 
angles are shown in the top left of Figure 2, again plotted with asterisks for SWEA and 
diamonds for STE. For SWEA, the blue perpendicular cut measurements are measured 
with bins with a pitch angle of 90°, while the green parallel cut measures electrons with 
pitch angle of 37°. We do not use the pitch angle closest to 0°, as these electrons will likely 
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lie inside the loss cone and, according to the fast Fermi theory, will not be accelerated by 
the shock. The lack of pitch angle coverage from the STE instrument during this time 
interval limits the ability to choose near-parallel and near-perpendicular cuts. This is 
mitigated by the fact that the quiet time superhalo is nearly isotropic [Lin, 1998]. 

The blue and green solid lines on the incident energy spectrum plot in Figure 2 are the 
results of fitting the eVDF to a model including the core, the halo, the strahl, and the 
superhalo. Each individual population is shown as a dotted line. The blue (perpendicular) 
solid line is the sum of the core, halo, and superhalo populations. The green (parallel) 
solid line is the sum of the blue line and the strahl population. 

Due to the previously mentioned instrumental charging effects, SWEA measurements 
of the core population are limited. To characterize the core, we instead use convected 
electron measurements from the SWE instrument on board the Wind spacecraft Ogilvie 
et al. [1995], calculating the time lag according to the method described in Opitz et al. 
[2009]. The time-shifted core density n c is 4.72/cm 3 , and the core electron temperature 
T c is 8.82 eV. This density corresponds to a local plasma frequency of ~ 19.6 kHz, which 
matches the frequency of the observed LFR Langmuir waves. 

The halo component for this interval is well respresented by an isotropic k distribution. 
The strahl component is modeled by subtracting the near-perpendicular SWEA cut from 
the near-parallel SWEA cut, then fitting the excess parallel observation as a function of 
energy using a spline fit. The superhalo is fitted with a modified K-type distribution which 
accounts for the relativistic behavior of the highest energy electrons [Xiao, 2006; Pierrard 
and Lazar, 2010]. It is assumed that the quiet-time superhalo is isotropic. 
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The model eVDF matches the data well over the range of the halo, strahl, and superhalo. 
The parameters for the fitted eVDFs are summarized in Table 1. 

4. Predicted foreshock electron properties 

As stated above, the aim of this paper is to compare observed foreshock electron obser- 
vations to the canonical Fast Fermi theory of electron acceleration. In order to use this 
theory to predict the properties of a foreshock electron beam, information about (a) the 
shock parameters, (b) the interaction of the electrons with the shock, and (c) evolution 
of the beam is necessary In this section, we discuss the estimation of these quantities. 

4.1. Shock parameters 

According to Fast Fermi theory, the acceleration of an electron by the bow shock depends 
on the loss cone angle a, which is determined by the magnetic mirror ratio between the 
upstream magnetic field and the field in the shock front, the HT speed wht, and the 
cross shock potential $. In order to determine the eVDF at an arbitrary point in the 
foreshock region, it is also necessary to know the depth of the spacecraft in the foreshock 
region, since the depth determines the cutoff electron speed v c necessary for an electron to 
reach the spacecraft (electrons with v < v c will be convected downstream by the motional 
electric field of the solar wind and will not be observed at the spacecraft [Anderson et ai, 
1979; Filbert and Kellogg, 1979; Cairns, 1987]). 

None of these parameters can be directly measured when the spacecraft is far from 
the shock. On 20 December 2006, STEREO/B was approximately 100 i?E from of the 
Earth's bow shock, and several days had passed since its previous bow shock crossing. 
Furthermore, the location of that previous crossing was not close to the quasiperpendicular 
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connection point at which the electron acceleration takes place. For these reasons, a and 
$ cannot be accurately estimated using data from the previous STEREO shock crossing. 
Previous papers dealing with electron acceleration at the bow shock (e.g., Filbert and 
Kellogg [1979]; Larson et al. [1996]) have used the straight field line approximation to 
determine 9b n , which yields ^ht, and depth in the foreshock, which yields v c . However, 
with a spacecraft more than several tens of Re upstream of the bow shock, magnetic field 
line wandering makes such an assumption invalid [Zimbardo and Veltri, 1996]. We must 
therefore estimate the relevant parameters using different methods. 

The ratio of the maximum field in the shock layer (including any overshoot) to the 
upstream field B maK /B u determines the opening angle a of the fast Fermi loss cone. We 
use data from the Geotail spacecraft to estimate this parameter. On 20 December 2006 at 
approximately 15:00 UT, Geotail observed several crossings of the quasiperpendicular bow 
shock at GSE coordinates of approximately [2.7, — 17.6, 3.8]i? E , i.e., in the neighborhood 
of the point where the STEREO/B spacecraft was connected to the shock some hours 
earlier. Although this conjunction is only approximate in both time and space, the Geotail 
data nevertheless give a picture of the typical conditions of the quasiperpendicular shock 
during the relevant time interval. Figure 3 shows magnetic field data from the Geotail 
MGF instrument [Kokubun et al, 1994] on 20 Dec 2006. The quasiperpendicular shock 
crossings are seen as rapid changes in \B\. The ratio B nmx / B u is measured for several 
shock crossings and shown above the plot. The ratios range from under 4 to slightly over 
5, with an average B max /B u of 4.64. Using the relation sin 2 a = B n / B max to compute the 
loss cone angle yields a = 27.7°. 
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In order to estimate $, we use the result of an analytical model [Kuncic et al, 2002] 
for the cross shock potential which depends only on the upstream perpendicular electron 
temperature and the magnetic compression ratio. 



This model has been shown to be consistent with spacecraft measurements of electron 
temperature made during bow shock crossings [Hull et al, 2000]. Using the convected 
Wind electron measurements, the Geotail magnetic field, and Equation (2), we estimate 
the cross shock potential $ m 72. V. 

The v ht parameter must be estimated from the foreshock electron data. Since electrons 
gain energy proportional to 2^ht, the energy gain of the electrons can be used to estimate 
wht, using the Liouville mapping technique described in the following section. We use the 
relatively low energy SWEA electron energy bins to estimate vht, by matching observed 
fluxes with predictions based on the upstream measurements. For the 07:25 UT interval 
at the foreshock edge, the estimated i>ht is 1500 km/s. For the 07:27 UT interval, the 
estimated Vht is 500 km/s. Similarly, v c can be estimated from the foreshock data by ob- 
serving the lowest energy bin which shows evidence of foreshock electrons. The estimated 
values of v c for the 07:25 foreshock edge interval and the 07:27 interval are 13000 km/s 
and 4250 km/s, respectively. 

The estimated Fast Fermi parameters a, $, vkt, and v c are listed in Table 2. 

4.2. Liouville mapping 




(2) 
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Liouville's theorem states that, in the absence of collisions, the phase space density / of 
a collection of particles remains constant along a trajectory defined by the (conservative 
and differentiable) forces which govern the system. In the case of the solar wind electrons, 
this implies that /; for the incident eVDF remains constant along trajectories determined 
by the fast Fermi theory, so long as that theory remains applicable. The theorem allows 
for the prediction of foreshock eVDFs given the upstream observations and the shock 
parameters estimated above. 

The Liouville theorem is the basis behind previous theoretical work predicting foreshock 
eVDFs [e.g. Cairns, 1987; Fitzenreiter et al, 1990]. It has also been used in experimental 
studies analyzing the evolution of the eVDF through bow shock crossings [Gosling et al., 
1989; Hull et al, 2001; Lefebvre et al, 2007]. In the HT frame, electrons traveling toward 
the shock which satisfy the condition 



will be reflected by the shock, assuming conservation of the adiabatic invariant \i (see, 
e.g., Knock et al. [2001, Eq. 9]). In the HT frame, these electrons reverse the sign of 
their parallel velocity, which adds a parallel velocity of 2-u H t- A schematic diagram of fast 
Fermi acceleration which combines the illustrations of Leroy and Mangeney [1984] and 
Wu [1984] is shown in Figure 1 of Pulupa et al. [2010]. 

In order to be observed at the spacecraft, the electrons must additionally meet the 
criterion that their parallel velocity is higher than the geometric cutoff velocity at the 
spacecraft. The Liouville mapping of the incident eVDF fa to the foreshock eVDF fa can 
then be summarized in the following equation: 




(3) 
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v\\ < max[w c ,f H T] 



ff(v\\,v±) = < 



fi(v\\,v ± ) v\\ > max[u c ,UH T ], 



(4) 



in loss cone 
fi(2vHT-V\\,v±) v\\ > max[u c ,VH T ], 

reflected 



Using Equations 3 and 4, we can directly calculate the predicted foreshock eVDFs given 
by Fast Fermi theory, and observe deviations from the measured STEREO data. 

4.3. Quasilinear relaxation 

Foreshock electron beams can create a bump on tail feature which is unstable to the 
growth of Langmuir waves on very short time scales. The foreshock beam is unstable 
when the reduced parallel distribution f\\(v\\) displays a positive slope in the v\\ direction 
(df\\/dv\\ > 0). In the quasilinear theory of beam relaxation [Grognard, 1975; Muschietti, 
1990; Robinson and Benz, 2000], positive slopes in the eVDF are mediated by Landau 
resonance, which generates Langmuir waves. In short order, the region around the bump 
on tail evolves to a marginally stable (df\\/dv\\ « 0) plateaued state. 

In this work, the method of Knock et al. [2001] and Kuncic et al. [2004] is used in order 
to estimate the region over which the flattening of the beam will occur. The procedure is 
as follows: Using the foreshock eVDF obtained from the Liouville mapping, the reduced 
distribution f\\ is constructed. The bump on tail is flattened by solving numerically for 
the points V- and v + which satisfy the conditions that f\\(v-) = f\\(v + ) and the number 
of particles JJ+ f\\dv\\ is conserved. The velocity of the beam v\> and its width Afb can 
then be calculated from i>_ and v + . 

5. Foreshock electrons 
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Energy spectra and PADs for the foreshock electrons are shown in the middle and 
bottom rows of Figure 2, with the foreshock edge observations in the middle row and the 
deeper foreshock observations in the bottom row. As in the top row, energy spectra are 
to the left, PADs are to the right, SWEA data is plotted with asterisks, and STE data 
with diamonds. The foreshock beam is evident in the foreshock PADs as an enhancement 
(compared to the solar wind PAD) for angles greater than 90°. Near-perpendicular and 
near-antiparallel (backstreaming) cuts through the eVDF are plotted in the left panels 
with blue and red solid lines. The solid lines plotted along with the electron spectra are 
the results of the Liouville mapping. 

We analyze two intervals in order to emphasize different aspects of the foreshock electron 
beam. Using the foreshock edge observations, we can isolate and examine the behavior 
of the high energy suprathermal foreshock beam electrons, lower energy electrons having 
been cut off from the beam by the foreshock geometry. Using the deeper foreshock ob- 
servations we can observe the lower energy foreshock beam. These energies are useful for 
analyzing the effect of the strahl population and the limitations of the fast Fermi theory. 

5.1. Foreshock edge 

The foreshock edge observations from 07:25 UT are shown in the middle row of Figure 
2. For reference, the incident /strahl population line from the top row is also plotted 
on this spectrum as a green solid line. A foreshock beam consisting of electrons which 
were simply reflected without energization would lie directly atop the green line. Electron 
measurements showing evidence of energization therefore lie above the green line. The 
electron beam in the foreshock edge energy spectrum is evident from the SWEA energy 
bin at 1.1 keV up to the 20 keV STE energy bin. In terms of the Liouville mapping, 
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this implies that the velocity cutoff is somewhere below 1.1 keV, and that i>ht can be 
determined by matching the observed beam at 1.1 keV to Liouville map predictions. For 
the foreshock edge, we find v c ~ 13000 km/s, and i>ht ~ 1500 km/s. However, it is clear 
that the Liouville predictions do not match the STE observations at higher energies. 

This result is consistent with observations of Gosling et al. [1989], who noted that a Li- 
ouville mapping predicting the downstream eVDF differed significantly from observations. 
Additional processes beyond /i-conserving reflection must apply for these suprathermal 
energies. 

Figure 4 shows only the parallel STE data from the foreshock edge energy spectrum of 
Figure 2. The solid red line is the predicted spectrum using the fast Fermi model, as is 
also shown in Figure 2. The red shaded area in Figure 4 brackets the model predictions 
with extreme values of vwv- Varying t>HT between and 10000 km/s yields the lower and 
upper boundaries of the shaded area, respectively. It is clear that both the level and shape 
of the observed spectrum differs from the model predictions for the low and high values 
of i>ht, emphasizing that the departure of the observations from fast Fermi is not simply 
a consequence of parameter choice. The data are fit to an power law with an exponential 
roll-off, following Oka et al. [2006]. For this particular interval near the foreshock edge, 
we find a relatively energetic power law, with a power law index Y of 2.04 and a roll-off 
energy Eq of 5.5 keV. The power law index and rollover energy evolve as a function of 
time with variation in the electron beam, and the relationship between T and E Q may be 
the subject of a future study. 

5.2. Mechanisms for high energy electrons 
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Previously, Gosling et al. [1989] noted that observed electron fluxes could not be ad- 
equately explained by the fast Fermi mechanism, and suggested that the suprathermal 
upstream electrons could be leaked from a downstream population which had been accel- 
erated by the shock. For the mechanism for generation of the observed enhancements in 
suprathermal electrons, Gosling et al. [1989] proposed shock drift acceleration. Shock drift 
theory has been shown to be equivalent to the fast Fermi mechanism [Krauss- Varban and 
Wu, 1989] for single encounters, but also can be naturally extended to study multiple en- 
counters between the electron and the shock [Ball and Melrose, 2001]. Shock drift theory 
[Vandas, 2001] has been used to study the parametric dependence of upstream electrons 
on shock parameters, and to compare observed and predicted spectra. The results are 
qualitatively similar but suggest that additional processes are important, especially for 
the high energy (above 10 keV) electrons. 

More recently, several alternative ideas have been proposed, which explain the enhanced 
suprathermal electrons as a consequence of some type of process which can create pitch 
angle scattering in an incident electron [Burgess, 2007]. The scattered trajectory en- 
ables the electron to make multiple encounters with the shock and be accelerated to high 
energies. We briefly describe several proposed processes in the paragraphs below. 

Recent simulation work has also considered the effect of non-ideal structure on shock 
fronts on suprathermal electron acceleration. Using hybrid simulations and an injection 
profile of 100 eV electrons, Burgess [2006] investigated several effects of shock front ripples, 
finding a considerable increase of the electron maximum energy and electron acceleration 
from a broader range of ^bn than was the case for planar shocks. The difference in the in- 
teraction between a rippled and a planar shock lies in the ability of electrons to encounter 
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the shock multiple times, due to the connection of a magnetic field line to multiple points 
on the surface of a rippled shock. Evidence of such multiply connected electron accel- 
eration sites (although for larger length scales than studied by Burgess [2006]) has been 
observed at CME-driven interplanetary shocks and is related to the generation of Type 
II radio bursts [Bale et al, 1999; Pulupa and Bale, 2008]. 

Previous Geotail observations of energetic electrons show a connection between ener- 
getic electrons and upstream whistler waves [Oka et al., 2006], showing a dependence of 
foreshock electron power law index on the whistler critical Mach number. This dependence 
has been explained as accelerated foreshock electron beams generating wave activity which 
can scatter electrons back and forth across the shock in a mechanism akin to the diffusive 
shock acceleration observed in the ion foreshock [Amano and Hoshino, 2010]. In a single 
event analysis using data from a Geotail shock crossing, a gradual (as opposed to spiky) 
temporal profile of accelerated electrons combined with observed upstream whistlers has 
been observed [Oka et al, 2009]. 

Guo and Giacalone [2010] investigated the effects of upstream turbulence, also using 
hybrid simulations. Using test particle trajectories, Guo and Giacalone [2010] directly 
demonstrated the effect of multiple shock crossings on electron energization. Guo and 
Giacalone [2010] also found that the level of upstream turbulence affected the maximum 
attainable energy for an injection profile of 100 eV electrons, as higher levels of turbulence 
increased field line wandering and thus the effect of multiple shock encounters. 

Savoini and Lembege [2001] and Savoini et al. [2010] have used two dimensional par- 
ticle simulation to directly simulate electron acceleration at the terrestrial bow shock. 
By tracing test particle orbits, several different types of accelerated electrons are found, 
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including a population of electrons consistent with Fast Fermi and populations which are 
temporarily trapped in the front or escaped downstream. The latter populations make 
up the higher-energy portion of the foreshock beam observed in the simulations, and pos- 
sibly correspond to the high energy power law tails observed in the hybrid codes and in 
observations. 

Recent work on quasiperpendicular shocks shows evidence of nonstationarity driven by 
reflected ions on the time scale of the ion gyroperiod causing the shock front to self- 
reform dynamically [Krasnoselskikh et al, 2002; Lobzin et al, 2007; Lembege et al, 2009]. 
Lembege and Savoini [2002] have suggested nonstationarity as an explanation for the 
bursty nature of upstream electron events. The nonstationarity on the shock front, which 
in simulation results generates a highly non-uniform shock surface, could provide a mech- 
anism for the multiple encounters which can generate high energy electrons. 

In future work, observations such as those presented in this paper may be of use in 
discriminating between the various mechanisms described above. With the STE detector, 
parameters of the suprathermal electron spectrum such as T and E can be studied with 
good energy resolution up to the highest beam energies. Systematic variation of these pa- 
rameters with solar wind and shock parameters may vary differently among the proposed 
mechanisms described in this section, and comparison with observations may therefore 
offer insight into the physics of electron acceleration. 

5.3. Limit of Fast Fermi 

The foreshock edge beam examined above exhibited a high velocity cutoff. As a con- 
sequence, the lower energy bins were not observed in the foreshock beam. Deeper in the 
foreshock, the foreshock electron beam extends to lower energies. For the observations 
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from 07:27 UT, shown in the bottom row of Figure 2, the electron beam is apparent 
in the energy spectrum from below 100 eV and up to the first several STE energies at 
several keV. Using the same method of estimating velocity parameters as used in the 
previous section, we find v c ps 4250 km/s, and i>ht ~ 500 km/s for the deeper foreshock 
measurements. 

Using the parameters from Table 2, we can again compare the fast Fermi predictions 
with the observed data. The Liouville-mapped data (red line) fits the observations (red 
asterisks) well for the relatively low energy SWEA bins from 151 to 400 eV, with an energy 
beam which is slightly enhanced compared to the incident strahl population from the top 
row (green line). This enhancement is consistent with fast Fermi acceleration. However, 
an additional enhancement in the foreshock electron beam beyond anything consistent 
with Fast Fermi theory begins in the 650 eV SWEA energy bin. 

For typical conditions at the terrestrial bow shock, the combined SWEA and STE 
energy range covers the transition between adiabatic and non-adiabatic electrons, as the 
electron gyroradii for electrons with energies of tens of eV are generally much smaller than 
the shock scale (defined by the ion inertial length or ion gyroradius [Bale et al, 2003])., 
while the electron gyroradii at tens of keV are comparable to or greater than these lengths. 
Using convected Wind measurements from the SWE Faraday cup instrument, it is possible 
to estimate the ion Larmor radius and the ion inertial length. For this event, ru = 74 km 
and c/cupi = 123 km. Setting the electron Larmor radius to be equivalent to these length 
scales requires electron energies of 87 keV and 32 keV, respectively. Similarly, an electron 
with an energy of 650 eV has a gyroradius equivalent to 0.086 r Li or 0.143 c/u pi . Since 
these ion scale lengths define the scale size of the shock, this confirms that the electron 



DRAFT 



February 17, 2012, 4:16am 



DRAFT 



PULUPA ET AL.: STEREO FORESHOCK ELECTRONS X - 23 

foreshock beam begins to diverge from Fast Fermi at the expected point: when is no 
longer much smaller than the scale size of the shock ramp and hence /i is not conserved. 

We note that Gosling et al. [1989] used a Liouville technique to map the upstream 
eVDF through the shock and compare to the observed downstream eVDF. In Figure lib 
of Gosling et al. [1989], the comparison between the observed and Liouville-mapped eVDFs 
shows the same deviation from single encounter fast Fermi theory, with the observed eVDF 
substantially enhanced beginning at energies of several hundred eV. 

5.4. Role of the strahl 

Since the cross shock potential at the terrestrial bow shock is typically of the order 
of the core electron temperature or higher [Hull et al, 2000, 2001; Kuncic et al, 2002; 
Lefebvre et al, 2007], the foreshock electron beam is constituted primarily of electrons 
from the suprathermal solar wind populations. Simulations of electron acceleration and 
the resulting plasma wave activity [Wu, 1984; Krauss-Varban et al, 1989; Krauss-Varban 
and Burgess, 1991; Kuncic et al, 2004; Malaspina et al, 2009] model the upstream eVDF 
as a sum of core and halo populations, where the halo population is represented as either 
a hotter Maxwellian population or as a k distribution. 

The effect of the geometrical velocity cutoff is to separate the foreshock eVDF into 
incident solar wind electrons v\\ < v c and reflected electrons v\\ > v c (See Equation 4). 
In the case of a backstreaming beam consisting of reflected halo electrons, a significant 
bump on tail will form only if the incident electrons are strongly accelerated by a high 
i>ht at the acceleration point — that is, only in a narrow region where #bn is very close to 
90°. 
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In the case of a backstreaming beam where the incident eVDF has a strong strahl 
component, however, the bump on tail can be more easily formed. If the ratio of strahl 
density to halo density n s /n h is high, even a small energization will result in a significant 
bump when combined with the effects of the geometric cutoff — the geometric cutoff will 
produce an eVDF which consists of the incident halo below v c but the reflected strahl 
above v c . A wider range of 9^ n connection points will be able to produce bump on tail 
eVDFs with a strong strahl component present. 

This effect is illustrated in Figure 5. Using the eVDF parameters from Table 1, the 
fast Fermi parameters for the deep foreshock from Table 2, and Equation 4, the reduced 
distribution f\\ computed for two cases: one where the strahl population is included in 
the incident eVDF, and one where it is removed. In each bump on the tail of /y 

is created by the fast Fermi process. In the case where the strahl beam is included in the 
eVDF, the bump is higher and the beam is wider. To quantify the difference in the bumps, 
the width of the beam At>b is calculated as described in Section 4.3. A rough estimate 
for the energy in the electron beam is also calculated, by comparing n^v^ normalized to 
n^v^. The beam from the eVDF including the strahl population is more than twice as 
wide in parallel v and has more than three times the beam energy as does the beam from 
the eVDF without the strahl (see Figure 5). 

The geometry of quasiperpendicular connection to the bow shock divides the electron 
foreshock into two separate wings, and that the strahl population is incident on the shock 
in only one wing. Given the discussion above, the presence of the reflected strahl in one 
foreshock wing will lead to an asymmetry in foreshock wave activity. Such asymmetries 
have been noted, and correctly explained as a result of strahl electrons, using statistical 
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studies of wave activity in the foreshock of Venus [Crawford et al, 1993, 1998]. Statistical 
analysis of electron and wave measurements has shown that the asymmetry also exists 
in the terrestrial foreshock [Pulupa et al, 2011]. Recent results also connect observations 
of Langmuir waves in solar wind magnetic holes to the presence of a strahl population 
[Briand et al, 2010]. Briand et al. [2010] note that the suprathermal nature of the strahl 
population allow it to enter the magnetic hole and provide a source population for electron 
beams. Similarly the suprathermal nature of the strahl places the strahl electrons above 
the threshold condition for reflection determined by the cross shock potential $. In both 
the foreshock and in magnetic holes, the strahl adds asymmetry to the eVDF. When the 
asymmetric eVDF interacts with a magnetic mirror-like structure such as a shock or the 
edge of a magnetic hole, the asymmetry leads to an enhancement in the bump on tail 
feature which is responsible for the generation of Langmuir waves. 

It is also possible that the microphysics of the generated Langmuir waves may be af- 
fected by the presence of the strahl. Recent observations of type Ill-related Langmuir 
waves [Malaspina et al, 2011] have pointed to a relationship between electron beam speed 
and wave polarization. By supplying a relatively high-energy seed population to the ac- 
celeration process, the strahl may indirectly affect wave polarization. 

6. Summary and Conclusions 

We have presented observations of foreshock electron beams using the STEREO space- 
craft, combining the SWEA and STE electron instruments to make observations which 
simultaneously cover the entire range of foreshock electrons with good energy resolution 
and sensitivity. Analysis of the beam electrons leads to the following conclusions. 
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Our results are consistent with the canonical Fast Fermi picture of foreshock electron 
acceleration for low energy electrons. A significant enhancement of electron acceleration 
begins at energies where the electron gyroradius is a significant fraction (~ 0.1) of the ion 
gyroradius and/or ion inertial length. The high energy tail of the foreshock electron beam 
is well fitted by a power law with an exponential roll-off, and extends up to several tens 
of keV before falling below the level of the background superhalo electron population. As 
discussed in Section 5.2, a number of possible mechanisms have been shown to produce 
suprathermal electrons in this energy range. 

Presence of the strahl is a significant factor in producing the bump on tail eVDFs which 
generate Langmuir waves in the electron foreshock. In particular, in upstream regions 
where the cutoff velocity v c lies in the halo-strahl energy range, the asymmetry generated 
by a strong strahl implies that a large Vrt is not necessary to produce a significant bump. 

The experimental results presented here may be useful for differentiating among the 
several physical mechanisms of the shock acceleration of electrons. We suggest that simu- 
lation and modeling work in the future should include incident electrons featuring realistic 
strahl and superhalo populations in addition to halo electrons. 
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STEREO Behind 




Figure 1. Upstream foreshock electron event observed by STEREO/Behind on 20 December 
2006. From top to bottom, the plotted quantities are: GSE magnetic field, electric field spec- 
trogram from S/WAVES LFR, energy spectrogram from the STE D3 detector, STE PADs at 
13.6 keV and 5.2 keV, and SWEA PADs at 1060 eV and 400 eV. The foreshock electrons are 
visible in the STE spectrogram starting at around 7:25 UT. The PADs for the electrons show 
enhancements in the antiparallel direction corresponding to foreshock electron beams. The time 
interval delineated by the two vertical bars near 7:20 UT is used to calculate the incident eVDF. 
The time intervals near 7:25 and 7:27 UT correspond to the two foreshock intervals analyzed in 
Section 5. 
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Figure 2. Energy spectra and PADs for the incident electrons (top row) and the foreshock 
electron beams at the foreshock edge (middle row) and deeper in the foreshock (bottom row). 
SWEA pitch angle bins are marked by asterisks, and STE bins are marked by diamonds. A full 
description of this figure can be found in Sections 3 and 5 of the text. 
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Figure 3. Several crossings of the quasiperpendicular bow shock as observed by Geotail on 
20 December 2006. Geotail crossed the shock near the electron acceleration site for the event 
observed by STEREO/B. The crossings are indicated by dotted lines. For each analyzed crossing, 
the ratio B m3jX /B u is shown. A full description of this figure can be found in Section 4.1. 
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Figure 4. STE energy spectrum for the foreshock edge electron beam, demonstrating significant 
enhancement of accelerated electrons compared to the Fast Fermi theory. The spectrum is well 
fit by a power law with an exponential roll off (black solid line). The red line and the red shaded 
area show the range of predicted spectra calculated using Liouville's theorem and the fast Fermi 
acceleration mechanism. A full description of this figure can be found in Section 5.1. 
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Figure 5. Reduced parallel distribution /n from the Liouville mapping. The solid line in the 
upper plot shows f\\ calculated using the parameters described in Sections 3 and 4.1. The large 
bump on tail is a result of the strong strahl population in the incident eVDF. The horizontal 
dotted line represents the relaxed plateau distribution constructed by eliminating the bump on 
tail while conserving particle number, as described in Section 4.3. The lower plot shows the 
results of the same procedure with the strahl removed from the incident eVDF, with the result 
that the bump is significantly smaller. A full description of this figure can be found in Section 
5.4. 
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n c 


4.72 cm" 3 


T 

± c 


8.82 eV 


n h 


0.065 cm" 3 


T h 


48.6 eV 




5.9 




0.731 xl0~ 6 cm 




2.50 keV 




2.7 



Table 1. Measured parameters for the incident eVDF. The field-aligned strahl component is 
fit empirically using a spline fit and not to a parameterized model. 



a 27.7° 
$ 72. eV 

v ht (foreshock edge) 1500 km/s 
vht (deeper foreshock) 500 km/s 

v c (foreshock edge) 13000 km/s 

v c (deeper foreshock) 4250 km/s 
Table 2. Parameters used to calculate the Liouville mapping of electrons through the fast 

Fermi shock encounter. These parameters are estimated using upstream plasma measurements 

and Geotail shock crossings. See Section 4.1 for a complete description of parameter estimation. 
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